Theory for absorption of ultrashort laser pulses by spheroidal metallic nanoparticles 
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The theory for the electric and magnetic fields energy absorption by small metallic particles 
subjected to the irradiation by ultrashort laser pulses of different duration in the region of surface 
plasmon excitation is developed. For the particles of the oblate or prolate spheroidal shape there 
has been found the dependence of the absorbed energy on a number of factors, including a particle 
radius, a degree of the shape deviation from a spherical one, a pulse duration, the orientation of 
the magnetic field upon particle, the magnitude of carrier frequency, and the value of a shift of the 
carrier frequency of a laser ray from the frequency of the surface plasmon excitation in a spherical 
particle. An appreciable absorption grows at the length of free electron pass large compared to 
the particle size is established. The phenomenological and kinetic approach is compared each with 
other. 

PACS numbers: 78.67.Bf; 68.49.Jk; 73.63.-b; 75.75.+a 



I. INTRODUCTION 



Metallic nanoparticles (MNs) are interesting objects 
from the viewpoints of physics of condensed state and 
practical applications. An increase in the local electric 
fields in the close vicinity of nanoparticles allows the uti- 
lization of the MNs as biological markers^ and can to lead 
to entirely new types of biological sensors^ The mod- 
ern optical investigations make it possible to test the 
optical response of a separate nanoparticle and, thus, 
to study the properties of separate nanoobjects. 3 This 
opens the new opportunities for the direct probing of 
the many-electron dynamics in confined systems. The 
remarkable results were obtained recently in works 4-£ 
In recent years, the ultrafast dynamics of electrons in 
MNs attracts the permanent attention of experimenters. 
Nanostructures in general have been widely used in mod- 
ern devices of high-speed electronics and optoelectronics. 
Owing to the recent achievements in the ultrafast spec- 
troscopy technique, it has become possible to monitor 
the femtosecond dynamics of the electron gas confined in 
separate MNsii The use of short-duration pulses makes 
it feasible to study the dynamics of fast processes occur- 
ring in atoms, molecules, and solids. Pico- and femtosec- 
ond resolutions allow to study of vibrational and rota- 
tional intramolecular motions, the dynamics of carriers 
in semiconductor nanostructures, the phase transitions 
in solids, the processes of formation and breakdown of 
chemical bonds, etc<&2 The possibility to generate pow- 
erful femtosecond pulses has lead to the appearance of 



a new, rapidly developing branch 
high-energy physics. 
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When a MN is placed in the field of a monochro- 
matic electromagnetic (EM) wave and the wavelength is 
far greater than the particle size, the absorption is ob- 
served associated both with the separate and collective 
kinds of electron motion. The input to the separate ab- 
sorption gives as an electric component of an EM wave 
(electric absorption), so the magnetic component (mag- 
netic absorption) as well. The collective electron mo- 
tion become apparent in so-called surface resonances con- 
nected with the excitation of plasmon vibrations (plas- 
mons) of an electron gas i 10 ' 11 There are three such reso- 
nances in an ellipsoidal particle and only one in a spher- 
ical particle. In the case of a monochromatic wave, the 
collective contribution becomes dominant provided that 
the wave frequency is close to one of the frequencies of 
plasma resonances^ Numerous applications in decora- 
tive painting^ and optical devices^ are based just on 
the resonances of surface plasmons. However, in spite of 
the appreciable theoretical efforts^ a body of problems 
has remained unsolved hitherto. 

In previous worksi^ we had shown, that by exposure 
of metallic nanoparticle to EM wave, the frequency of 
which is far from plasmon resonances, a separate mecha- 
nism becomes actually in absorption. Depending on the 
shape and size of the particle, as well as on the frequency 
and polarization of wave, it was shown that, the mag- 
netic absorption can be higher or less than the electric 
one. The situation is drastically change with using of 
an ultrashort laser pulses for irradiation, because they 
contains almost all harmonics. The shorter the pulse, 
the greater is set of its harmonics. Thus, the ultrashort 
pulses can excite all the plasma resonances. This is as- 
sociated with the fact that, in such an ultrashort pulse, 



one can always find the harmonics which correspond to 
the resonance frequencies of surface plasmons. From the 
other side, the frequency distribution is very nonuniform 
(Gaussian) upon the spectra of an ultrashort pulses. How 
these factors can change the relative input of an electric 
and magnetic components in absorption, we intend to 
study below. Another aspect influencing the absorption 
is associated with the shape of a particle. Earlier, we 
showed that, in the case of nonspherical MNs, the optical 
conductivity, which determines both the electric absorp- 
tion and the half-width of plasma resonances, becomes 
the tensor quantity^ As far as we know, such a feature 
in absorption of ultrashort laser pulses by nonspherical 
MNs was not investigated in details. 

In this paper we study for the first time, the pecu- 
liar features of the absorption by the MNs of ultrashort 
pulses, deduced by an electric and magnetic components 
of EM wave. Such a problems remains less studied, es- 
pecially for the particles of a nonspherical shape, whose 
dimensions are less or exceed the electron free path in it, 
and have attracted a lot of interest of experimenters in 
recent years. In particular, we will find: (i) the energy 
absorbed by MNs from ultrashort laser pulses of a differ- 
ent duration; (ii) its dependence on the particle shape, 
frequency of a carrier wave, pulse duration, and parti- 
cle sizes; (iii) the shift of the resonance frequency of the 
absorption with respect to the plasmon frequency in a 
spherical MN. 

The rest of the paper is organized as follows. The 
model and an initial foundations of problem are presented 
in Sec. II. Section III contains the study of a plasmon res- 
onances in the MNs. In section IV the study of magnetic 
absorption is demonstrated using both the phenomeno- 
logical and the kinetic approaches, and Sec. V contains 
the conclusions. 



II. MODEL AND STARTING POSITIONS 

Consider a case where an MN is irradiated by a laser 
pulse, whose electric field is given as 



E(r,i) = E exp 
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cos 



u>o\t- 
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wo 



(1) 



where T is the quantity reciprocal to the pulse duration, 
cjo is the carrier frequency of an EM wave, |ko| = ojq/c, 
and Eo is the maximal value of an electric field in a 
pulse. In addition to the electric component, the field 
of a laser pulse contains also the magnetic one which 
is associated with the former through the corresponding 
Maxwell equation 



rot E(r,£) 
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c at 
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We set E(r, t) in the form ([T} and will find the magnetic 
component of a pulse, by using Eq. @ . Relation ^ takes 



the simplest form if one operates with the Fourier com- 
ponents of the above quantities. Actually, for an electric 
field, e.g., we obtain from Eq. $Q 
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Carrying out the Fourier transformation in Eq. <j2j) , one 

finds connection between Fourier components of the mag- 
netic and electric fields as 



H(r,w) 



ILC 



-rot E(r, w). 



Using the last part of Eq. ([3]) , we find for magnetic field 



H(r,w) = ^(mxEo) 



x I e «^ + e JF 2 
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where m = ko/fco is the unit vector directed along the 
EM wave propagation. 

The Fourier component of the electric field amplitude 
takes the simplest form in the limit T — > 0: 

E(r, w) r ^ - TrEoe* 01 "^ [% - w ) + S(u + wo)]- (5) 

Whereas the electric field of a laser pulse creates the 
potential electric field Ei n (r, t) inside an MN, the mag- 
netic field induces the eddy current field E c d(r, t). To 
calculate the energy absorbed by an MN, we should know 
the inner fields E in (r, t) and E c d(r, t). To find these 
quantities, we pay attention to the following feature of 
the coordinate dependence of the Fourier components 
(|3]) and (|4|): when the characteristic dimension R of a 
nanoparticle is such that 



fc -R« 1, 



(6) 



i.e. the wavelength of the carrier wave is far greater 
than the MN dimension, the coordinate dependence of 
the Fourier components E(r, u) and H(r, w) can be ne- 
glected inside the particle. This means that, in order to 
determine the inner fields inside an MN, provided that 
inequality (JSJ) is satisfied, we can set the Fourier compo- 
nents of these fields such as they would be in the spatially 
uniform fields: 

E(r,w)->-E(0,w) and H(r,w) ->■ H(0,w). (7) 

In the case of an nonspherical MN which has, for exam- 
ple, the form of an ellipsoid, this allows us to write the 
components of the inner electric field as 
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and the magnetic eddy current field as 
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by analogy with the procedure described in Ref. [16]. 
The other components of the eddy current field can be 
obtained from Eq. (|9]) by means of a cyclic permutation 
of indices. The notations in formulas jSJ) and (j9|) are as 
follows: Lj (j = x, y, z) are the geometrical factors, also 
known as the depolarization factors^ e(cj) = e'(oj) + 
is"(uj) is the frequency-dependent complex permittivity 
of the particle, R x , R y , R z are the ellipsoid semiaxes in 
the x, y, and z directions, respectively. 

The linear dependence of the eddy current field on co- 
ordinates can be easy seize from Eq. ©, which deter- 
mined it. If one perform the time Fourier-transform of 
Eq. @ , than the right-hand side of this equation [as it 
is seen from both the Eq. ((4]) and condition ©], one can 
consider as permanent, coordinates-independent quan- 
tity. This means that rot E ec ;(r,u;) = const. The last 
equality is satisfy only if 



[E ed (r,u)]j = ^2ctjk(u)xk, 



(10) 



fc=i 



e.g., the curly field depends on coordinates linearly. Here 
a is some matrix does not depended on the r; its com- 
ponents will be further specified. 

It is seen from expressions © and © that the uni- 
form external electric field induces also the uniform elec- 
tric potential field inside the ellipsoidal MN, whereas the 
uniform external magnetic field generates the coordinate- 
dependent eddy current field. The inner fields Ei n (r,t) 
and E ec ;(r, t) induce the currents inside the MN with the 
corresponding densities jm(r, t) and j ed (r, t). As a result, 
the particle absorbs the energy of the EM field of an in- 
cident laser wave. The total absorbed energy is sum of 
electric and magnetic absorption and can be presented 

„„17 
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x [Rej e (r,i)ReE* n (r,i) 

+ Rej cd (r,i)RcE * d (r,i)], (n) 

where the integration should be carried out over the 
whole particle volume. Here, W(i) is the absorbed power. 
Depending on whether the absorption is brought about 
by E; n (r, t) or E ed (r, t), it is called the electric or mag- 
netic absorption, respectively. In view of the equality 
j*(w) =j(-w) and 
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which follows, in particular, from the Parseval relation for 
the Fourier integral^ expression (|11[) can be rewritten 
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x [j e (r,w) E* n (r,w) +j cd (r,w) E* d (r, u) + cc.](12) 

The values of Ej n (w) and E c d(r,w) can be calculated 
from formulas ([SJ and (|9]), respectively. Thus, the prob- 
lem remaining to be done is to find the Fourier compo- 
nents of the current densities j; n (r, oj) and j e d(r, w). In 
the general case, the current produced by the inner fields 
E; n (o;) and E d(r,w) at a point r of the particle can be 
expressed as the integral over all values of electron veloc- 
ities v 



j(r,w) =2e 



/ m \ : 



d 3 vv /i(r,v,w). 



(13) 



Here, v — |v|, /i(r, v, w) is the Fourier component of 
the nonequilibrium distribution function usually consid- 
ered as an addition to the equilibrium Fermi distribution 
function /o(e) which depends only on the electron ki- 
netic energy e. Usually, the way to calculate the function 
/i(r,v, ui) consists in the solution of the corresponding 
linearized Boltzmann kinetic equation. As a rule, this 
equation is written for the timc-dcpcndcnt distribution 
function (see, for example, Ref. Ilq ). Performing the 
Fourier transformation of this equation with the use of 
expression Q and 



/i(r,v,w) 



we obtain equation for f\\ 



/i(r,v,t) e iut dt, 



(v-iu)fi(r,v,u) +v— /i(r,v,w) 

+ ev[E in (w)+E ed (r,u;)] |-/ (e) = 0, (14) 

where v is the collision frequency in the particle bulk. To 
move further, we should add the corresponding boundary 
conditions to Eq. (TT4"|) . We take assumption of the diffuse 
reflection of electrons from the inner surface of a particle 



/i(r,v,w)| s =0, 



<0 



(15) 



as the boundary conditions. Here, v n is the velocity com- 
ponent normal to the surface S. The substantiation of 
such a boundary conditions and the solution of (|14l) are 
presented, in particular, in Ref. Ila . 

Below we consider separately the peculiarities of the 
electric absorption when the plasmon resonances are ex- 
cited by ultrashort laser pulses and the magnetic absorp- 
tion take place due to magnetic component of the laser 
wave in more details. 



III. PLASMON RESONANCES 

To this point, we have considered the general approach 
which includes both the electric and magnetic absorp- 
tions. It was stressed above that the absorption can be 
either separate or collective. When an MN is in the field 
of a EM wave and the wave frequency is far from the 
plasmon resonance, the separate mechanism of absorp- 
tion prevails. In this case, as was shown in Ref. Ilia . 
the electric or magnetic absorption can be dominant de- 
pending on a number of factors, namely, on the wave fre- 
quency, its polarization, the MN size and shape. In the 
present section, we concentrate, however, on the features 
of the absorption of ultrashort laser pulses. In this case, 
the main contribution to the absorption results from the 
plasma resonances, i.e. from the collective mechanism, 
which is a constituent part of an electric absorption. For 
this reason, we will take into consideration in present 
section only that part of absorption in Eq. dill) , caused 
by the electric component of the EM field. To find the 
total energy of electric absorption, we will start from ex- 
pressions (|T2j) and (|13[) . and determine the form of the 
function /i(r, v,i), by solving Eq. (|Ti)) with the corre- 
sponding boundary conditions. Then, we will substitute 
the obtained function into Eq. (fT5|) and calculate the 
Fourier components of the current density. Finally, by 
taking into account that the Fourier components of the 
inner field are described by expression ([5]), we will sub- 
stitute the obtained Fourier components of the current 
density into Eq. (|12[) with E e d = 0. The similar proce- 
dure was described in detail in Ref. Qjl and, thus, we 
won't repeat it. We only note that the kinetic equation 
(fT4f almost completely coincides with Eq. (27) of work 15 
with the only difference: Ei n (w) should be changed by 
Ei oc (o;). What is more, the integrand of the integral over 
frequencies in Eq. (|12l) formally coincides with the initial 
expression in Ref. U]& Thus, we can use the results of 
Ref. Ila . in particular formula (75), and write the total 
energy of an electric absorption in the form 
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where the frequency dependence of E, (0, w) is given by 
expression ©. Here, 
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is the half- width of a plasmon resonance, which strongly 
depends on the geometrical shape of a particle and is 
defined in terms of diagonal components of the optical 
conductivity tensor a. For the different directions j, 
it is governed by the aforementioned factor Lj. For a 
spheroidal particle (R x = R y = R±), these components, 
for example, are reciprocal to the square of the frequency 
at high frequencies and have the form 
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where vf is the electron velocity on the Fermi sphere, 
uj p is the plasma frequency of electron oscillations in a 
metal, and T](e s ) and p(e s ) are some smooth functions 
depending on the spheroid eccentricity e s . The analytical 
dependence of these functions on the ratio between the 
spheroid semiaxes as welll as, e s can be found in Ref. Ila , 
For a spherical particle, r\ = p = 2/3. 

The half-widths of plasmon resonances are an im- 
portant characteristic, since it contains the informa- 
tion about the character of interactions in the system. 
Through the interaction between surface plasmons and 
single-particle excitations, the electron scattering gives 
rise to the decay of surface plasmons, which reveals itself 
in the broadening of lines of the differential transmission 
spectra. The surface plasmon decay can occur in different 
ways, depending on the particle size. Often, the contri- 
bution of the Landau decaj*^ to the linewidth, which is 
dominant in the case of a small particle radius a of the 
order of 5-20 A^2£ is calculated. At greater particle sizes, 
the Landau decay competes with the decays caused by 
other interactions. Recently, we studied 2 ^ how the rela- 
tion between the transverse j±_ and longitudinal 711 com- 
ponents of the plasmon resonance half- width depends on 
the degree of ellipsoid's oblateness or prolateness for the 
frequencies that are higher or lower than the characteris- 
tic frequency of electron reflection from the particle walls. 

In the case where the plasmon decay is not significant, 
we can replace the Lorentzian in expression (|16[) by a 
(^-function as^ 2 - 
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Then, with accounting of Eq.Q, the expression for the 
absorbed energy takes the form 
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where 



(19) 



(20) 



We limit the analysis of Eq. (|T9j) to the case of a 
spheroidal MN and choose, for convenience, such a po- 
larization of the electric field, which makes it possible to 
excite the plasmon oscillations of electrons both along 
and across the spheroid rotation axis. If the field com- 
ponents Eg II and Eoi_ are directed along and across the 
rotation axis, respectively, then such a polarization is, in 
the general case, 
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FIG. 1. Dependence of the plasmon resonance frequencies on the 
form of gold nanoparticle. 

where 9 is the angle between the spheroid rotation axis 
and the vector of the incident electric field strength E . 
Then, frequencies ([20)1 can be rewritten as the frequen- 
cies of electron oscillations along and across the axis of 
rotation: 



L 



ll.-L 



(22) 



which, for the given R\\ ± will be the frequencies of the 
corresponding plasmon resonances. The dependence of 
the plasmon resonance frequencies on the particle form 
is illustrated in Fig. 1 for Au. 

We will use the electromagnetic radiation energy w 
traversing the nanoparticle during the total duration of 
the pulse as a scale. This energy at normal light inci- 
dence, we define as 



E o (0,t)| 2 dt = — J-(l + e ^ 



where c is the light velocity, and 



Eo| 2 , 
(23) 

(24) 



is the particle section. In what follows, for the sake of 
illustration, we will calculate the quantity 



2w 



(25) 



which is the ratio of the energy absorbed by a unit of 
the MN volume to the energy traversed the nanoparticle. 
The ratio between particle volume and particle section is 
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where 
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FIG. 2. Dependence of the energy absorbed by a spheroidal gold 
particle with R = (R^ijii) 1 ' 3 = 200A at the plasmon resonance 
frequency on the ratio between the spheroid semiaxes for different 
values of V, s" 1 : 2.637 X 10 15 (1), 1.582 X 10 15 (2), 1.13 X 10 15 
(3), and 0.879 x 10 15 (4). 8 = 7r/4. 

Figure 2 shows the dependence of S e on the degree of 
ellipsoid's oblateness or prolateness in the case where = 
7r/4 and the frequency equals that of surface plasmons of 
a spherical particle uj = uj p /V3 = tt. Here and below, 
the calculations are carried out with the use of formulas 
(IT91 and (|2~T1) - (|2"7)l . given electron density in the Au 



particle A n sa 5.9x10 
1.37 x 10 16 s" 1 andfi 



22 cm~ 3 



j p = (Anne /m) 



1/2 



7.91 x 10 5 s . The numerical 



data for Y are taken here and below (in illustrations) 
from the Table 1. It is convenient to measure T in the 
fractions of the carrier frequency wo- Then, the duration 
of a laser pulse can be set in the carrier wave periods 
or wavelengths, that is l/r„ = n/ujQ ~ n\o, where n 
is an arbitrary integer. Thus, for the region of plasmon 
resonances (when loq = Q), each wavelength, which is a 
multiple of A U7 can be brought into correspondence with 
a certain pulse duration. 

TABLE I. Correspondence between the value inverse to 
the pulse duration and the wavelength of a carrier wave 
in the region of plasmon resonances for Au particle in 
vacuum. 



r x io 1 



A 



7.911 
2.637 
1.582 
1.130 
0.879 
0.791 
0.527 
0.396 



Ao 

3A 
5A 
7A 
9A 
10A 
15A 
20A 



As follows from curves 1~4, the laser pulses of longer 
duration are absorbed more effectively. As for the shape 
dependence of the absorption, the maximum of absorp- 
tion is achieved at R±/R\\ = 1, i.e., when the particle is 
in vacuum and have the spherical shape. If the particle is 
embedded in the dielectric matrix with e > 1, then max- 
imum is shifted to the red side of the spectra range, pro- 



portionally to the value of e. For example, for Au particle 
embedded in the glass matrix with e — 7, the plasmon 
resonance pick occurs at the frequency ll m ~ 11 — 0.5511. 
The calculation of S e was also carried out by means of 
the numerical integration in formula (|16[) for the various 
particle radii R within the range of R = (RuR^) 1 ^ 3 — 
50 -7- 400 A. On the whole, it turned out that the value 
of S e , calculated numerically, is (i) slightly smaller than 
the values obtained on the basis of formula (IT91 and 
(ii) depends on R very weakly for various 24 T. At the 
same time, the latter dependence is more pronounced at 
smaller T provided that R is not great (R = 50 -j- 150 
A). For example, for particle with R — 100 A at 
T = 0.879 x 10 15 s -1 , the highest possible values of 
S e calculated analytically with the use of Eq. (flUj) ex- 
ceed those numerically determined from formula (|16p by 
about 7%. For greater R and T, this percentage de- 
creases. For this reason, we will rest in the further calcu- 
lations of S e on the evaluations performed according to 
formula 03 
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FIG. 3. The same as in Fig. 1, for the frequency luq = 5.4 X 
10 15 B -1 . 

As soon as the carrier frequency deviates from that 
of the surface plasmon in a spherical particle, the peak 
of absorbed energy depending on the ratio R±/R\\ splits 
into the two ones: the first maximum corresponds to the 
prolate particles (Rx/R\\ < 1), whereas the second - to 
the oblate particles (R±/R\\ > 1). This fact is illustrated 
in Fig. 3 for the frequency oj — 5.4 x 10 15 s _1 and 
various values of the incident pulse duration. 

If the pulse duration is fixed, then, changing the fre- 
quency of a carrier wave, we can trace the evolution of 
the dependence of the maximum of absorbed energy on 
the MN shape. Figures 4(a) and 4(6) illustrate such 
a dependencies. For comparison, the absorption curve 
at a frequency Wo = fi, where H is the resonance fre- 
quency of a spherical MN, is shown (curve 1 in each 
part of the figure). For the frequencies exceeding H 
(see Fig. 4(a)), the absorption maximum first shifts to- 
wards the side of more oblate particles with increase in 
wo and then splits into the two peaks: the first one is 
observed at R±/Ru > 1, whereas the second one - at 
R±/R» < 1. Figure 4(6) demonstrates the data calcu- 



lated for the frequencies lower than 11, with the calcula- 
tions carried out for the same values of the pulse dura- 
tion. As the frequency gets lower, the absorption max- 
imum shifts towards the side of more prolate particles 
and, afterwards, splits into two, as in the former case. 
As ujo deviates from the resonance frequency, the peak 
of the absorbed energy first decreases in absolute value, 
then splits and finally, as the carrier frequency strongly 
differs from 11, the values of both emerged peaks reach 
the plateau. 
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FIG. 4. The energy absorbed by a spheroidal gold particle with R = 
(-R^-Rn) 1 ' 3 = 200A as the function of the degree of its oblateness 



or prolateness for T 



0.879 X 10 1 



at the frequencies ojq 



higher than that of the plasmon resonance $1 of a spherical MN 
(a): u) = 8.8 X 10 15 s" 1 (2), 9.7 X 10 15 s" 1 (3), and 10.6 X 10 15 
s" 1 ( 4) and lower than Q (b): u) = 7.0 X 10 15 s^ 1 (2), 6.1 X 10 15 
s" 1 (3), and 5.2 X 10 15 s" 1 (4). 9 = tt/4. Curves 1 in both parts 
of the figure correspond to the case luq = Q. 

These results are in compliance with the correlation, es- 
tablished by us earlier,— between the geometrical shape 
of a nanoparticle and the frequency of the light ab- 
sorption, depending on them, the rapid growth in the 
strength of light pressure on a particle is observed. 

Let us choose the prolateness or oblateness of a 
nanoparticle fixed and trace how the absorption changes 
as a function of the carrier frequency. In the case of 
a prolate MN, with increase in the duration of the in- 
cident laser pulse, one observes not only the growth of 
the peak of absorbed energy, but also its splitting (as 
r becomes higher than a certain critical value) into two 
peaks (a doublet) located on the opposite sides of the line 



ujq/SI = 1 (Fig. 5), where fl corresponds to the resonance 
frequency of a spherical particle. It is worth to note that 
the peaks at ujq < tt and ujq > il correspond to the plas- 
mon resonance frequencies wm and wj_, respectively (see 
Eq. (J22I) and Fig. 1), whereas the inverse situation is 
characteristic of oblate particles: the former peak corre- 
sponds to lu± and the latter - to wii (see Fig. 1). 
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FIG. 5. Dependence of the energy absorbed by a prolate (Rj_/Rn = 
0.7) gold particle with R = (R^R^) 1 / 3 = 200A on the frequency 
of a carrier wave for different values of T, s — 1 : 1.582 X 10 15 (1), 
0.791 X 10 15 (2), 0.527 X 10 15 (3), and 0.396 X 10 15 (4) (taken from 
Table 1). d = tt/4. 

Thus, contrary to the case of spherical particles, the 
spheroidal MNs, either prolate or oblate, exhibit a dou- 
blet in the spectrum of absorption provided that T is 
sufficiently small. The doublet originates from two plas- 
mon resonances excited in particles by ultrashort laser 
pulses. At 6 = 7r/4, both the peaks of the doublet 
have the same height, which, being practically indepen- 
dent of a degree of the particle prolateness or oblateness, 
strongly depends on the pulse duration. In the case of 
prolate MNs, a minimum of the dip between the absorp- 
tion peaks is achieved at frequencies uiq lower than f2 
and shifts towards the longer waves with increase in the 
prolateness degree. The other tendency is characteristic 
of oblate MNs: the minimum is observed at frequencies 
higher than $7 and, as the oblateness degree increases, it 
first shifts towards the shorter waves and then, at suf- 
ficiently high oblateness degrees, the movement changes 
its direction. 

The peaks of absorption display the similar behavior. 
As the degree of the MN prolateness increases, the dis- 
tance between the peaks of the doublet grows due to 
the more drastic shift of the peak, caused by the plas- 
mon resonance at the frequency wii , toward the longer 
waves. At the same time, the height of the peaks re- 
mains constant provided that the degree of prolateness is 
small enough. Quite similar behavior is also character- 
istic of oblate MNs: the distance between the peaks of 
the doublet grows with increase in the degree of the MN 
oblateness. The reason for this is either the more dras- 
tic shift of the peak, caused by the plasmon resonance 
at wii, toward the highfrequency side of the spectrum 



(at sufficiently small degrees of oblateness) or the reverse 
movement of the peak, caused by the plasmon resonance 
at uj_ (at higher degrees of oblateness). The height of 
the peaks remains constant as the degree of oblateness 
keeps small enough. 

Under various experimental conditions, a small shift, 
either blue or red, with respect to the value u>o = tt 
has been observed. A number of microscopic approaches 
has been developed to explain such a shift (see, for ex- 
ample, Ref. 111] ). The most successful approaches are 
thought to be those built on the basis of a jelly model 
with the use of the theory of linear response in the lim- 
its of time-dependent local density^ However, these ap- 
proaches don't account for the shift caused by a change in 
the particle shape. It is the most important feature of our 
approach that results in the facts that the shape of the 
MNs specifies the frequency, at which the resonance ab- 
sorption occurs, and the shift of the plasmon absorption 
peak is associated with a change in the particle shape. 
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6. The energy absorbed by a prolate (Rj_/Rn = 0.7) gold 
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particle with R = (R^R^) 1 / 3 

quency of a carrier wave for V 



200A as the function of the fre- 



0.396 X 10 



ir, 



and different 



angles of the laser pulse incidence: 8 = 
dotted line, and 9 = n/6 - dashed line. 



7r/4 - solid line, 8 = n/3 



The relative height of the peaks in the doublet can be 
tuned by means of changing the angle of the laser pulse 
incidence. Figure 6 allows one to trace how the plasmon 
absorption along and across the spheroid rotation axis 
changes as the angle of incidence 9 is reduced from 7r/3 
to 7r/6. For convenience, we consider the prolate particle 
with i?j_/i?|| = 0.7, as in Fig. 5. As is seen in Fig. 6, the 
peak intensities are the same at 6 = 7r/4, but they be- 
come different as 9 deviates from this value. So, for exam- 
ple, at 9 — 7r/3, the peak associated with the resonance 
at uju is characterized by a reduced intensity, contrary to 
the peak associated with wj., which is characterized by 
an increased intensity. At 9 — 7r/6, however, the inverse 
situation is observed. It is obviously that, whereas the 
intensity of the peak at wii reaches a maximum at 9 = 
and turns into the zero at 9 = 7r/2, the peak at uj± has 
maximum intensity at 9 — ir/2 and minimum one - at 
6> = 0. 

Lastly, we consider the qualitatively different limit 
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transition when the width of an incident laser pulse is 
sufficiently large. To do this, we use the next represen- 
tation for (/•-function 2 ^ 



6(x) 



— = lim — exp(— x la ), 
Jit a->o a 



and calculate, according to Eq. ([3]), the value of 
\Ef\0,uj)\ 2 , which is included in Eq. (JTBJl 



,(o) 



2E o 3 



'7T\3/2 1 



[2) f;[S(u - uo) + S{u + uo)]. 

(28) 
Now relations (jTSJ and (28) yield 



^r(o,^) 



2F /7T\3/2 1 

f e |r-o = — (7?) fL 



Wo«7-jj(w )-Eo- 



(29) 



As is seen from Eq. (|29[). the absorbed energy comparing 
with Eq. ([TO1) is proportional only to the first order of 
the pulse duration 1/r. Expression ([29]) divided by the 
pulse duration gives us the value of the energy which is 
absorbed by a nanoparticle on average in a unit of time 
at its irradiation by long-duration pulses. This quantity 
differs by the factor 2 y7r/2 from the analogous quantity 
obtained for the absorption of a spherical wave^ 



The energy absorbed by spheroidal MP in the case of 
broad pulses is given from Eqs. (|2T)]) . (fT8"|) by 
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w e \r^o 



9 \V \ 2 V 
32 J^v 



2, ,2 v f 



"o^p 



V(e s 



cos" 



p(e s )sin 2 6> 



(lo 2 - c 2 ) 2 + [2 7 ||K)] 2 ^ H - ojI) 2 + [2 7± ( Wo )] 2 w 2 , 



(30) 



r 



where 



7( j)( w o) = 27rL/|N(T/|x(cjo), 



(31) 



and L11, L± are the geometrical factors along and across 
the long spheroid axis of a revolution, correspondingly. 



IV. MAGNETIC ABSORPTION 

If the frequency of a monochromatic EM wave is far 
from a plasmon resonance, then a separate mechanism 
plays the main role in the MN absorption. As we have 
already pointed out in this case, either an electric or mag- 
netic absorption can dominate. But if the size of a spher- 
ical particle, for instance, is larger than 50 A, then the 
magnetic absorption begin to exceed the electric one even 
at frequencies less than the frequency of electron reflec- 
tions from opposite walls of a particle?^ In this section, 
we will consider in detail, the behavior of the magnetic 
absorption at frequencies close to the plasmon resonances 
of the MN. 

In accordance with Eqs. (TTTJ) and (|12[) . the magnetic 
field energy can be presented in an integral form as fol- 
lows 



1 |>00 7 P 

Wm= o/ o~ / dr[j eri (r,u;)E* d (r,w) + c.c.]. (32) 
2j-oo 2tt J v 

The value of E ec ;(r,w) is already known for us from 
Eq. (5). In order to calculate j e d(r, w) with the use of 



the formula (|13|) , it is necessary first to find the distribu- 
tion function for electrons. For this purpose, we would 
seek for solutions of Eq. (JTJ|, where E in (uj) = 0. To sat- 
isfy both the equation (fT4|) and the boundary conditions 
(fT5j). we should pass to the deformed variables 



R 
Hi 



R 
Rj 



Then an ellipsoid can be represented as a sphere of the 
same volume with the radius R, and the solution of Eq. 
(|14p we can write in the form 

/i(r',v',w) = -^-^-^aijMv'jRjRi 

Q \ l _ e -(v—kv)t'(r'y) 



R 2 de 



' d(y — iui) 



V — UjJ 



where 



i'(r',v') = -=[rV + ^{R 2 - r' 2 )v' 2 + (r'v' 



(33) 



(34) 



In Eq. (|33|) . the energy derivative of the distribution func- 
tion in zero approximation in a small ratio of ksT /ep can 
be replaced by 



dfo 
de 



-6(s-sf), 



(35) 



where ep it the Fermi energy. The diagonal components R y = R±, R z = R\\). Using Eqs. ([9]) and (l38l) . the sum 



of the matrix o.jj are 0^=0. The nondiagonal ones can 
be expressed in terms of appropriate components of the 
magnetic field. For example, 



it. 



Rl 



Rl 



HP{0,u). 



(36) 



The other two components of the matrix a can be ob- 
tained from Eq. (1361) by means of a cyclic index permu- 
tation. The remaining three components one can find 
using the antisymmetric character of the a, namely, tak- 
ing into account the property a xy (ui) = —a yx (uj). 

Thus, to find the energy of the magnetic absorption 
it is necessary to undertake step-by-step the following 
procedure: to inset the obtained function /i(r, v, t) into 
Eq. (|13[) and then, taking into account the expression 
(HJ, to put the Fourier components for both the current 
density and the curly field into Eq. (p?2"j) . 

One can escape the above procedure for energy calcu- 
lation, using the phenomenological approach for an es- 
timation of current. It means that the absorption of a 
particle can be written down in a general form with real 
and imaginary parts of a dielectric permeability, and thus 
it is not necessity to solve the kinetic equation (fT4"]) with 
the boundary conditions (TTSl) . 



A. Phenomenological approach 

Assume that the particle size is much more greater 
than the electron free pass within it. Then, the curly cur- 
rent can be expressed in terms of a curly field by means 
of 

jed(w) = CTmM EedM, 

and the expression (|32[) can be rewritten in the form 

Wm = oY, / KM + ^mM] o- / |E J ed (r,w)| 2 dr. 

(37) 
One can change here the integration over an ellipsoid vol- 
ume V by the integration over the sphere of an equivalent 
volume. Then, it is easy to show, that 
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x * dr = V !±£ : 
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Rl 



j/ d - V T- I 
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z 2 dr = V^, 
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'V 

xy dr — I yz dr — I zx dv = 0. (38) 
v Jv Jv 

Further, we will restrict ourselves with consideration of 
the metal nanoparticles having the spheroid shape (R x = 



of integrals over a nanoparticle volume will be 



tj^^-jgy 



R 2 



H?,(0 jW )- 



tit |Jt_|_ 



R fl 



m 



Hi(0,w) 



(39) 



Here, Hi,(0,u) is the intensity of the magnetic field 
along the spheroid axis of a revolution, and H^(0,w) - 
transverse to it. It is worth to pay attention that the 
magnetic absorption depends on the magnetic field po- 
larization in the same way as the electric absorption de- 
pends on the electric field polarization. 17 Moreover, one 
can generalize the expression (|39|) to the case of an ar- 
bitrary coordinate systems if the field components are 
represented as 



H^ 



H^ H; 



H 



(H ■ n) n, 



(40) 



where n is a unit vector directed along the spheroid axis 
of a revolution. Then with an account of Eq. (pi!?)) , the 
expression (1371) transforms into 



[a m (u) + a* m (u)} (^) 
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(41) 



The field H(0,w) in Eq. (|4"Tj) is defined by the expression 
Q, and a m for a metallic nanoparticle in the phenomeno- 
logical case is 



M 



~7, 



4ir v 



iui 



(42) 



where uj p = Ait n e 2 /m is the plasma frequency of the 
electron vibrations in a metal. 

Next, we will perform calculations for certain polariza- 
tions of the magnetic component of the EM wave with 
respect to the particle orientation. We will consider the 
following two cases: 

(i) Let us direct the vector of magnetic field trans- 
verse to the spheroid axis. Then, the second item under 
the integral in Eq. (|4"T|) vanishes, and after substituting 
Eqs. (@| and (|4"2j) into Eq. (|4"Tj) . we get for the absorbed 
energy the expression 
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(43) 
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where 



f{u) = exp 



■wo) 



21 



AT 2 



exp 



(uj + UJ f 

4r 2 



(44) 

Our attempt to fulfil the integration in Eq. (j43|) in an 
analytical form exactly was not successful, but using the 
fact that the frequency of electron collisions with phonons 
v is a small quantity comparing to the frequencies uj s=s 
wo, which yield the main contribution to the integral, 
one can neglect by v under the integral in Eq. (|43|) . For 
instance, oj ^10 15 s _1 and for Au is °c ~3.39xl0 13 s _1 , 
that comes from the formula for electroconductivity of 
metal a = 1/p ~ ne 2 /(mi>), with ,oo°c = 2.04x 10~ 6 Q cm 
(Ref. [271 ) - Then, one gets the result equal to twice the 
integral of 



/(w) duj. 



(45) 



This integral can be calculated analytically even in the 
pointed limits. However, it is easy to show, that the 
result of the integration practically is not changed if one 
brings the lower integration limit in Eq. (|45[) down to the 
zero. This means only that the integral between limits 
[0, v\ 



|H(0,oj)| and substituting Eq. §4§ into Eq. (|4Tj) . we come 
to the next expression 



«w« = ~ (^r) |Eo| 2 /°° [aH + ^Mi/H^ 2 g. 

(49) 
Using both the Eq. (|4"2l and the result of the integration 
given by Eq. (|47p. we obtain finally the next expression 
for the energy absorbed with this polarization 



w„ 



1 V V 

8oT/2¥r 



u p R\ 
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I En 



(50) 



provided that u> 3> v. The square of an amplitude of the 
electrical field in the previous expressions and further for- 
mally could be changed by |E | 2 — > |H | 2 , because the 
maximal value of the magnetic field in the pulse would be 
specified like similar one in Eq. (TTJ), with taking into ac- 
count that the relation between fields is given by Eq. ([2} . 
If one suppose that a nanoparticle is irradiated by a 
plane EM wave, i.e., by the wave ([1]) with r —> 0, then the 
expressions (|48|) and (|50|) coincide asymptotically with 
the known results from our earlier work a 15 ' 16 , with an 
accuracy up to the constant 2y7r/2 [due to the distinc- 
tion in a pulse forms (see previous note 2 ^)]. 



B. Kinetic approach 



/(w) dtu < I v 



(46) 



is much less than one of I v , and one can neglect by it. 
Thus, for the integral between limits [0, oo], we have: 



exp 
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AT 2 

■■ V2tt T ( 1 + exp 
and as a result 
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exp 
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4r 2 



,2 1 



2r 2 



du> 



(47) 




(48) 



provided that u> 3> v. 

(ii) Now we pass to the case of EM wave-polarization, 
where the magnetic field is directed along the main 
spheroid axis. In this case the scalar product in a sec- 
ond item under the integral in Eq. (|4"Tj) is H(0,oj)n = 



When the nanoparticle size is less than the electron 
free pass in it, the phenomenological approach cannot 
be used anymore and ought to be changed by the ki- 
netic approach. Note, that the latter permits to obtain 
the correct results for the case when the particle size is 
larger than the electron free pass as well. In the kinetic 
approach, the current j e d(w) enters into Eq. (|3^|) . should 
be calculated by means of the formula (Tr"3"|) . Using the 
nonequilibrium distribution function (f3"3"]). one finds 
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' k x\ (l - e -("-^)*VX)) dr ' ; 



(51) 



where vf is the electron velocity on the Fermi sphere. We 
have omitted the item in Eq. (f5Tj) associated with the 
second term (with v[) under the sum sign in Eq. (|33|) . 
It is not difficult to show that the contribution of this 
term with integration over all electron coordinates will 
be equal to zero, because it is even with respect to co- 
ordinates, whereas an eddy field is an odd coordinate 
function. Then, the integration over all coordinates in 
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Eq. (|51j) . associated with the first term under the sum 
sign in Eq. ([3"3")l , can be done exactly. We would write 
here only the final result. The details of calculations can 
be found in Ref. [la. 



one of which will exceed, and another one will be less 
than the frequency associated with electron oscillations 
between particle walls 
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(52) 
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v a =v F /{2R). 



(55) 



The factor of 2 in Eq. (|54[) appear due to the fact that the 
integrants S(w) should be the even functions of uj. The 
case with uj < i>s will be named further as the lowfrc- 
quency (LF) case, and the case with uj > v$ - as the 
highfrequency (HF) one. Let us consider first the latter. 
As a basis for it, we will apply the approximations 
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(53) 



The next integration over all velocities space fails to be 
carried out in an analytical form for a common case. It 
can be fulfilled only if the frequency interval [-oo, oo] to 
split arbitrarily into the two parts 



oo us oo 

/ E(uj)dco =2 E LF (uj)duj + 2 E HF (io)dio, (54) 
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(56) 



tp 2 {v',uj) = il>a(v',w) - -ipi(v',u). 
Details of calculations are presented in Refs. [1^ and [l6|. 
Here, we write down only the final result which can be 
obtained after integration over all electron velocities with 
the account of Eqs. fl3D) and ([56]) 
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H.n.,;),,^^. (57) 
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In expression ([57| . pH{e s ) and »7m(e s ) are the functions 
depended only on the spheroid eccentricity e s . Their an- 
alytical forms are given in Ref. [T^. From Eq. ([5"7| . one 
can see that the frequency dependence of the integrand 
(for any given polarization) in the HF-case is contained 
only in the magnetic field amplitude. Practically this 
means that we can use the result of the integration ob- 
tained above in the phenomenological approach. If the 
lower integration limit in the phenomenological case was 
the value of z/, then vg will play here the same role. It 
is not difficult to ensure that in the case with z/j, one 
can diminish the lower limit of integration to the zero as 



well. Hence, the energy of a magnetic component of a 
laser EM wave in (i) polarization (when the vector of the 
magnetic component is directed transverse to the long 
spheroid axis of revolution) , with an account of the ex- 
pressions (01 and (|4"?| in Eq. ([57]) . finally becomes 
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(58) 



and for the (ii) polarization, when, vice verse, the vector 
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of magnetic field component is directed along this axis, 
is represented by 



9 V 
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32^/2^ 16 c 2 T 



1 + e 3?* p ff (e s )i?j.|E | 



(59) 

Eqs. (|58|) and ([59)) coincide each with other for a spher- 
ical particle, because r]^(0) = 4/^(0). Besides, if one 
takes into consideration that the power of the absorbed 
energy is W m = w m T, where w m is given by Eqs. ([55)1 
or (|59|). and the parameter T tends to the zero in W m , 
then these equations will transform (with an accuracy of 
the constant of 2a/7t/2) into the known ones from our 



previous calculation ;) 15 ' 16 for a plane wave. 

It remains to consider the LF-case. The calculation of 
Eq. (j5"Tj) in this case can be done with the use of the next 
approximation 
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(60) 



Furthermore, not complicated but cumbersome calcula- 
tions lead us to the next result 
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(61) 



where pL(e s ) and ^(e s ) are some smooth functions of 
the spheroid eccentricity e s . The behavior of these func- 
tions depending on the shape of the MN, which is speci- 
fied by the ratio of R±/Ri i, one can find in RefKLa. For 
a spherical particle 



It can be done analytically. It is easy to ensure, how- 
ever, that similar to the calculation of the integral I in 
Eq. pp| . one obtains the same result using the integrant 
from Eq. (|63[) for II . In other words, one can get that 



p L (0)=p H (0)=2/3, ^(0) = 7^(0) = 8/3. (62) 



For an estimation of Eq. (fBTj) . it is necessary to calculate 
the integral 
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(63) 



W m S = W mLF + W mH F 



9 V 



w V F 



6AJ2^8c 2 r 



Rj 



h<^Iv 



(64) 



This means that in the sum of Eq. (py±|) , one can neglect 
by the integral of II ■ In this case, for an arbitrary angle 9 
between direction of the magnetic field and the spheroid 
axis of revolution, we will obtain finally the expression 
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(65) 



The results (|48|) and (|50|) obtained in the phenomeno- 
logical approach for a spherical nanoparticle go over into 
the ones described by Eqs. (|S"8|) and (|59"1) . or Eq. (|6"5"|) . ob- 
tained in the kinetic approach, if one performs formally 
the next replacement: 
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(66) 
with the ordinary condition that ui 3> v. In the case 
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of a spherical particle, the above substitution, obviously, 
takes the form: 



15 v f • i , , 

v — > , provided that 

16 R ' 



to » v. 



(67) 



The calculation are performed assuming the next 
numerical parameters for Au particle: the electron 
concentration 23 of 5.9xl0 22 cm" 3 , ?L = 7.91 x 10 15 s _1 , 
v = 3.39 x 10 13 s _1 , and v F = 1.39 x 10 8 sm/s. 



C. Results and Discussion 

Let us illustrate graphically the expressions obtained 
above analytically. We will calculate the ratio 



*~>m — 
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(68) 



between the energy absorbed by a unit volume of the 
MN and the energy traversed the nanoparticle, given by 
Eq. ([23]). 

Studying the dependence of optical properties of 
nanoparticles on their shape, it would be worthwhile to 
compare the absorption for particles of various shapes 
but with equal volumes. The transformation of a particle 
shape can be described by the ratio of radiuses R±/R\ i. 
The condition of fixed particle volume (V — ^-R^R^ 



const) with given ratio of R±/R\ 
of R± or of i?| | . For instance, 



defines the quantity 
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(69) 



where R is the radius of a sphere of an equivalent volume. 
Let us study at first the dependence of S m on the par- 
ticle shape. In Figs. 7 and 8 for various polarizations of 
EM wave, the dependence of S m on the degree of oblate- 
ness or prolateness of a spheroidal particle is depicted (at 
the frequency of the surface plasmon cjo = w p /\/3 = fl, 
which was a resonant one for a spherical particle in the 
case of the electric absorption). The calculations for _L 
and || polarizations were carried out using the next for- 
mulae 



&m± — 



IOtt 



-JT - 



R\ 3 {Rx/R\\) 2/3 



1 + {RJR\\ 



&m\\ — 



vD? (R 

20^ 



(R±/R\\ 



,2/3 



(70) 



(71) 



for phenomenological case, and 
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FIG. 7. The energy absorbed by Au particle for polarization of 
a magnetic field transverse (||) to it axis of revolution vs degree 
of spheroid oblateness or prolateness. Curves are obtained at the 
frequency of the plasmon resonance uiq = £1, for different R, A: 200 
(curve 1), 250 (2), 300 (3). The phenomenological dependence for 
R=200 Ais given by triangles. 



The curves 1 — 3 correspond to the different particle 
radiuses, whereas the curve of phenomenological depen- 
dencies builded on the base of the formulae (|70|) and (|7Tj) 
is given by triangles for R = 200 A. Comparing curves 
1 — 3, it is clearly seen that the energy of laser pulses 
will be absorbed more strongly in Au particles of a larger 
volume. When the magnetic field polarization is directed 
along the MN rotation axis (|| polarization), the absorp- 
tion grows with enhancing of the particle oblateness (Fig. 
7), while for the transverse polarization of the magnetic 
field (_L polarization), the absorption reaches the max- 
imum for some values of the ratio of R±/R\ \ (Fig. 8). 
In the phenomenological case, this maximum is reached 
for particles of a prolate shape when R±/R\\ = l/\/2, 
and in the kinetic one - when Rj_/R^ | = 1/y/Tl. If one 
compare the magnitude of the absorption maximum for 
the same pulse duration at _L polarization, which has 
been obtained in the phenomenological and the kinetic 
approaches from formulae (j48|) and (f58|) . then one can 
get the relation 



S„, = T^v F n 2 ^p H (e s )(R ± /R ll ) 1 / 3 , (73) 
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FIG. 8. The same, as in Fig. 7 for _L polarization of magnetic field. 

From Eq. ([7^]) it follows that this relation increases 
with the reduction of a particle radius and does not de- 
pend on laser beam parameters. The same one can see 
from Eqs. (l70l) - ([73|) . where the ratio between absorbed 
and traversed energy in the case of magnetic absorption 
does not depend on parameter T. For Au particle with 
the radius, for instance, 200 A one can see (comparing 
curves 1 and those labeled by triangles in Fig. 8) that 
the absorption in the kinetic case is approximately of two 
times of magnitude higher. 

With the extension of a MN prolateness, the magni- 
tude of the absorption is decreased in the case of _L po- 
larization of the magnetic field, and increased - in the 
case of || polarization. The reason for this is that the 
number of electron closed orbits is increased with grow- 
ing of a prolateness of a particle in 1 1 polarization, and 
decreased in _L polarization of the magnetic field. 

Contrary to the electric absorption, in the dependence 
of the absorbed energy on R±/R\ i, there does not occur 
neither the shift of maximum (for _L polarization) nor 
an appearance of any others maxima with the deviation 
of a carrier frequency from the frequency of the surface 
plasmon Q. 

For the case when loq < $7, the intensity of mag- 
netic absorption slightly increases with decrease of loq 
and asymptotically approaches the constant value, when 
Lj >n [see Eqs. ([58]), <j59j), gg]), and ([50])]. This is true 
for both polarizations as well as for both the phenomeno- 
logical and the kinetic approaches. But depending on the 
duration of an incident pulse, a certain peculiarities ap- 
pears. Namely, the absorption intensity is increased with 
an enhancement of the pulse duration at different carrier 
frequencies for both examined here polarizations. This 
grows is more rapid at frequencies luq <C ft, but has a re- 
striction connected with the evident fact that the space 
size of a pulse c/T cannot exceed the length of the carrier 
wave Ao in vacuum: T/luq > l/(2ir) sa 0.16. 

In contrast to spherical particles, the magnetic absorp- 
tion in nonspherical ones drastically depends on magnetic 
field polarization. So, the absorption of a spheroidal MN 



increases for moderate prolate and decreases for oblate 
particles (comparing to the particles of a spherical shape 
of the same volume) at _L field polarization. It is the least 
for particles of an oblate shape. In the case of || polar- 
ization, vice versa: the absorption noticeably increases 
for an oblate and decreases for a prolate MN compar- 
ing to the spherical one. Similar trends holds true for 
phenomenological approach in both polarization as well. 

If one choose a nanoparticle of a spheroidal shape with 
some fixed ratio of R±/R\ | and will change only the mag- 
nitude of the carrier frequency, then one can establish 
using Eqs. <j58j), ([59]), (gH]), and ([50]) for different polar- 
izations of an incident field the following, (i) For MN of 
a prolate shape (R±/R\ | < 1) with an enhancement of a 
laser pulse duration at frequencies u>q <C fi, the maximum 
of the energy absorption increases and does not split into 
the two peaks as we obtain above for prolate Au particle 
in an electric fields with T exceeding some limit quantity 
(see Figs. 5 and 6). (ii) If the ratio of ujq/VL is small, the 
rise of magnetic absorption reaches the higher point for 
the least frequency ujq = vs- (hi) Comparing intensity of 
absorption for pulses of equal duration, but different po- 
larization, one can ensure that more rapid enhancement 
of the energy absorption by MN with increasing of carrier 
frequency occurs at the _L polarization of the magnetic 
field. 

For oblate particles, vice verse, an enhancement of the 
energy absorption by MN with decreasing of a carrier 
wave proceeds faster for the 1 1 -polarization of the mag- 
netic field, then for it _L polarization. The pulses of larger 
duration (with small value of T) are intensively absorbed 
at lower frequencies of a carrier wave in both polariza- 
tions. 

Thus, the spheroidal MN of an oblate or prolate shape 
comparing to the spherical one can absorb the energy of 
the magnetic field from ultrashort laser pulses more or 
less intensively depending on the magnetic field orienta- 
tion with respect to the spheroid axis of revolution. 

It would be interesting also to monitor the absorption 
by a nanoparticle of a fixed shape with volume chang- 
ing. This can be easy fulfilled for fixed ratio of R±/R^ 
with changing the radius of an equivalent sphere R in the 
previously obtained formulae. Such one dependence for 
nonspherical particles can be strongly differ from an anal- 
ogous one in the case of spherical particles and requires 
a special investigation. There is only problem, when the 
volume of particle is changed with simultaneous preserva- 
tion of its shape, the size of a particle in same directions 
becomes greater than the electron free pass in those di- 
rections. Similar study for an electric absorption we have 
performed in Ref. l2lL 

From expressions (gU and dSU, or (gl]) and ([53]) one 
can see only that the energy of the magnetic absorp- 
tion with the use of the phenomenological description 
increases quadratically with R±, whereas at the kinetic 
description, the growth proceeds linearly with the R±. 
Generally, the increase of magnetic absorption with R± is 
due to the enhancement of the number of electron closed 
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orbits with growing of a particle oblateness in both po- 
larization of the magnetic field. The different behavior 
of the two approaches with R± can be understood from 
above Eq. (|66p . where phenomenological parameter v is 
replaced by v F / R±_ in a kinetic approach. 

By making use of the above results, namely, Eq. (jT9j) 
for the electric absorption and Eqs. (1551) and ([5"9")) valid 
for the magnetic absorption in a high-frequency case, one 
can compare relative contributions for both polarizations. 



close to or move away from the fixed carrier frequency. 
That is important for the amount of energy absorbed 
by particle. Such an approach is very demonstrable and 
productive for a theoretical studding of the optical prop- 
erties of a MN. But for direct comparison of a theoretical 
predictions with an experiments it is necessary to have 
the data for a metallic nanoparticles of equal volumes but 
different in their shapes. 
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For spherical particle, using Eq. (1B"2"|) . one can get that 
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An estimation of the ratio of absorptions [from Eq. ([76])] 
for particle with the radius of R = 50A, for which the 
input of the magnetic field overcomes the electric one 
for a separate absorption, 21 shows that at the frequency 
ujq = £1 the absorption associated with an electric plas- 
mon excitation is considerably higher than the magnetic 
one. The same one can see, comparing the magnitudes of 
absorption shown in Figs. 7 and 8 for magnetic absorp- 
tion, with electric one (Fig. 2) for Au particle at plas- 
mon resonances. However, the ratio between magnetic 
and electric fields input can enhanced with the pulse du- 
ration, and for plane monochromatic wave will be equal 
to 1 or even exceed 1, depending on the particle sizes. 
The substantial rice of this ratio at plasmon resonances 
can be realized only in || polarization of magnetic field 
for highly oblate particles. For instance, for Au particle 
of R — 200A, we can reach the same input of magnetic 
component as for electric one (with T = 2.637 x 10 15 s -1 ) 
at R±/ Ru « 800. The input of a magnetic component in 
the absorption becomes comparable and even can exceed 
the electrical one, as well, if a carrier frequency loq devi- 
ates from the plasmon frequency Q, especially at small 
values of T. One can find from Eq. ([TBI that, e.g., for 



Au particle with R 

-115 



200A, and T = 1.13 x 10 la s -1 , at 
£1 - cj | = 5 x 10 15 s _1 , w m /w e ~ 1. 

In general, the absorption of energy from ultrashort 
laser pulses by spheroidal MN depends both on the par- 
ticle and the pulse characteristics. If the parameters of 
laser pulses are given, the absorption by MN depends on 
its volume (for fixed shape) as well as on the its shape 
(for fixed volume). We have especially interested in the 
dependence on the particle shape, because with changing 
of it, we can variate the position of plasmon resonances 
on the frequency scale. In other words, with changing 
the ratio of R±/R\\, the plasmon resonances may came 



V. CONCLUSIONS 

We develop the theory of the absorption of ultrashort 
laser pulses of a different duration and carrier frequen- 
cies by small metallic particles of a spheroidal shape. The 
cases when an electron free pass is greater or much less 
than a particle size are considered. There has been ana- 
lyzed the dependence of the absorbed energy on the de- 
gree of a particle oblateness or prolateness at the fre- 
quency of a plasmon resonance and at frequencies that 
are greater or less than it. The simple analytic expres- 
sions are obtained, which make it possible to study the 
plasmon and magnetic absorption of laser pulses depend- 
ing on the duration of these pulses, the magnitude of the 
carrier frequency, the particle shape, and the field polar- 
izations. The analysis of the dependence of the energy 
absorbed by a spheroidal MN on the degree of its pro- 
lateness or oblateness was carried out for different values 
of pulse duration and frequencies that are higher, lower, 
or equal to the plasmon resonance frequency. 

At the frequency of a carrier wave, which is equal with 
that of a surface plasmon, the maximum absorption was 
observed for spherical MNs. As soon as the carrier fre- 
quency deviates from that of the surface plasmon in a 
spherical particle, the two maxima appear in the depen- 
dence of the absorbed energy on the ratio of spheroid 
semiaxes: one of them corresponds to the prolate par- 
ticles, while the other - to oblate particles. As the fre- 
quency deviates from the resonance one, the peak of the 
absorbed energy at first decreases in an absolute value, 
then splits and finally stabilizes for particles being more 
and more prolate or oblate. 

For the case of spheroidal MNs subjected to the laser 
irradiation, with increase in the pulse duration, the peak 
in the dependence of the absorption on the frequency of a 
carrier wave also splits into two peaks located on opposite 
sides of the resonance (for a spherical particle) frequency. 
The double peaks originate from the resonances excited 
by ultrashort laser pulses at the frequencies of plasmon 
oscillations along or transverse to the spheroid rotation 
axis. A distance from the minimum in a well between the 
doublet peaks to the frequency of a plasmon resonance 
for a spherical particle provides the information about 
the degree of prolateness or oblateness of a nanoparticle. 

As the degree of prolateness or oblateness increases, 
the distance between the doublet components grows. At 
the fixed incident angle of a laser pulse, the peak height 
keeps constant, with its value being dependent only on 
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the pulse duration. It was considered how the change 
in the incident angle of a laser pulse directed along or 
transverse to the spheroid rotation axis effects the value 
of the relative absorption intensity. 

We have found that the absorption in nonspherical par- 
ticles drastically depends on the magnetic field polariza- 
tion with respect to the particle rotation axis. Namely, 
for MN of a spheroidal shape, with the direction of the 
magnetic field transverse to the spheroid axis of revolu- 
tion (_L polarization), the absorption is increased for a 
prolate and decreased for an oblate MN compared to a 
particle of a spherical shape of the same volume. The 
case of the magnetic field orientation along the axis of 
revolution of a particle (|| polarization) proved to be the 
contrary one: the absorption is increased for a oblate MN 
and decreased for a prolate one. 

It appears that the energy of laser pulses with the 
larger value of 1/r is absorbed by MN better than with 
the smaller one. In the case with increasing of a par- 
ticle oblateness, the magnetic absorption is increased 
for longitudinal (||), and decreased - for a transverse 
(_L) polarization of the magnetic field. Comparing to 
the MN of spherical shape, the magnetic absorption for 
nanoparticles of a prolate shape is increased with ris- 
ing of particle prolateness at _L polarization of magnetic 
field, reaching maximum for the ratio of spheroid hal- 



faxes R±/R\ \ ss 7/15, and diminishes for the larger pro- 
lateness. At 1 1 polarization, increasing in a prolateness 
leads to decrease in the absorption. 

In the kinetic approach, the maximum value of the 
magnetic absorption for a fixed pulse duration and at 
_L polarization of magnetic field can exceed by two times 
or more (for particles with R < 200A) the absorption, 
obtained in the phenomenological approach [at T = 0°C]. 

With reduction in a carrier frequency, the intensity of 
the magnetic absorption by MN of a prolate shape is 
increased for _L polarization of the magnetic field more 
distinctly than at it || polarization. As soon as it took 
place, the pulses with longer duration are absorbed bet- 
ter at the lowest frequencies of the carrier wave. For 
oblate particles, on the contrary, the faster enhancement 
of the energy absorption by MN with increasing carrier 
frequency occurs for 1 1 polarization of the magnetic field, 
than at it _L polarization. 

For frequencies of carrier wave close to the plasmon 
resonances in spheroidal MN, the electric absorption is 
considerably higher than the magnetic one at _L polariza- 
tion of magnetic field and can be comparable with one - 
at || polarization for highly oblate MN. 

A comparison between the results obtained in both 
the phenomenological and kinetic approaches have been 
advanced. 
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